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I. INTRODUCTION 

Atoms and atomic ions have long held a place at the 
very heart of precision measurement and metrology. For 
example, the second is defined in terms of the Cs hy- 
perfine structure interval, the Rydberg constant is mea- 
sured by spectroscopy of atomic hydrogen and the ratio 
of electron to proton mass is known from the oscillation 
frequencies of trapped atomic ions. The importance of 
atoms in this field lies partly in the detailed understand- 
ing that we have of atomic structure but also in the tech- 
nical capabilities that exist for preparing and manipulat- 
ing atoms and ions. In recent years, both the compu- 
tational methods for understanding molecules more fully 
and the experimental methods for producing and con- 
trolling them have advanced enormously. This has led 
to a surge of interest in using molecules for precision 
measurements, especially where they offer new proper- 
ties that are not available from atoms and atomic ions. 
For example, the rotational, vibrational and electronic 
structures within a molecule offer a wider range of co- 
existing frequencies than one finds in atomic systems. 
Moreover, polar diatomic molecules have a built-in cylin- 
drical symmetry, whilst more complex molecules can have 
a handedness - structural conformations that atoms can- 
not offer. In this chapter, we discuss some applications of 
molecules to current problems in precision measurement 
and we outline recent technical advances that make some 
of these applications possible. 



II. TESTING INVARIANCE PRINCIPLES 
A. Do fundamental constants vary in time? 

Recent measurements have shown that the expansion 
of the universe is accelerating, requiring the Einstein 
equations of cosmology to have a "dark energy" term, 
which was previously assumed to be zero[T]. This sur- 
prise, together with the ongoing search for a quantum 
theory of gravity, has led theorists to question some of 
the most basic assumptions in our physical model of the 
universe, including the assumption that the constants of 
nature are indeed constant over time. 

In atomic and molecular physics, two of these con- 
stants are particularly important. They are the fine 
structure constant a, and the electron-to-proton mass 
ratio jjL — rrie/mp. The Rydberg energy Ry sets the gross 
scale of electronic binding energy. Relative to this, the 
fine structure splittings are characteristically smaller by 



the factor and the hyperfine structure is smaller again 
by a further factor of order /j,, since /j, relates the mag- 
netic moment of the nucleus to that of the electron. Con- 
sequently, it is possible to search for a variation of a by 
comparing fine and gross structure at two different times. 
Similarly the variation of /i can be deduced from a fur- 
ther comparison of the hyperfine structure with one of 
the other energy scales. Molecules bring an important 
new dimension to this search [2^ by adding two more en- 
ergy scales: vibrational energy of order Ry^Jra^/M and 
rotational energy of order Ry{me/M), where M is the re- 
duced nuclear mass. Uzan f3] has recently reviewed both 
the theoretical framework and the experimental tests for 
variable constants. 

So far, the most productive experimental method has 
been to study astronomical spectra, which permit mea- 
surements of Aa/a and A/i//i, typically with a precision 
of 1 part in 10^, over enormous time intervals of order 
10 Gyr, giving uncertainties of order 10~^^/year in the 
average value of a/ a or At this level, some obser- 

vations seem to hint at a variation^ ,5j, while others do 
not [3 E]. In order to interpret the astronomical data 
one needs to know the sensitivity of each transition to a 
and /i, which requires input both from laboratory data 
and from numerical modelling. Using H2 in this way, 
Reinhold et al. [S] have found an average variation over 
12 Gyr of /i//Lt = (1.7 ± 0.5) x lO'^Vyear, which dif- 
fers from zero by 3.4(T. The ground state A-doublet of 
OH has also been observed at large redshift. The tran- 
sition frequency for this interval relative to its hyper- 
fine structure yields information on both a and /i [7]. 
When compared with the structure of OH today, these 
astrophysical measurements looking back 6.5 Gyr indi- 
cate with 2a confidence that < 2.1 x 10~^^/year on 
average and is consistent with zero. A laboratory exper- 
iment using a molecular decelerator to make slow OH |5] 
(there is more on deceleration in Sec.|V]) has recently im- 
proved our knowledge of these transition frequencies in 
the current epoch. This will allow further improvement 
in the astronomical data to yield an even better determi- 
nation of /i and a in the near future. Another frequency 
with high sensitivity to /i is the inversion splitting of NH3 
due to tunnelling. Analysis of quasar absorption spectra, 
comparing the inversion line to rotational transitions, has 
placed a limit fi/jj, = (-1 ± 3) x IQ-^^yr"! HU] on the 
variation of /x. This leaves the astronomical measure- 
ments in an intriguing but uncertain state. Are there 
systematic errors yet to be uncovered? Do the constants 
vary or not, and if they do, is the variation irregular in 
time or perhaps non-uniform in space? 

There is a good prospect that laboratory measure- 
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ments can help to clarify the situation. Although the 
time intervals available for comparison are short - years 
rather than giga- years - this method provides a promising 
alternative because recently developed frequency comb 
methods can link optical frequencies directly to the Cs 
standard, giving an absolute frequency accuracy below 
one part in 10^^ Even higher accuracy, approaching 
one part in 10^^ is becoming available in measuring rel- 
ative frequencies [T2]. Several molecular experiments of 
this type are now under way. It has been suggested [13] 
that the cold, trapped molecular ions or HD+ could 
be used to measure changes in ^ to a precision of 1 part 
in 10^^. A second laboratory-based proposal [Uj is to 
compare the inversion splitting measured in a slow, cold 
fountain of ND3 to an atomic reference. A third project 
is in progress to measure vibrational transitions in SFg 



B. Testing fundamental symmetries 

The electromagnetic forces that bind atoms and 
molecules together obey Maxwell's equations and the 
Dirac equation, as synthesised in quantum electrody- 
namics. This field theory has three important symme- 
try properties: it is invariant under space inversion (par- 
ity, P) charge conjugation (exchange of particles and an- 
tiparticles, C) and time reversal (T). These symmetries 
have profound experimental consequences. For example, 
the eigenstates of atoms and molecules have definite par- 
ity (unless there are degenerate conformations), leading 
to selection rules for radiative transitions. For similar 
reasons, atoms and molecules cannot have a permanent 
electric dipole moment (EDM) (barring degenerate con- 
formations). For example, the well known EDM of am- 
monia is not permanent because the nitrogen atom tun- 
nels back and forth at a frequency set by the splitting of 
the two opposite-parity field-free states. Of course, this 
splitting is small, and therefore it takes only a modest 
electric field to induce a dipole moment. 

For some years it was surmised that all interactions 
possess these symmetries, but an experiment in 1956 
[TB] showed that weak interactions violate parity, as seen 
by the fact that radioactive decay particles are emitted 
with a large left-right asymmetry. Within a decade, an 
experiment on the decay of kaon particles showed that 
strong interactions are also not symmetrical, having a 
small asymmetry under the combined operation CP [T7] . 
There is a theorem for the type of theories used to de- 
scribe particle interactions (local, Lorenz-invariant field 
theories) that they must be invariant under the triple 
reflection CPT. Since CP symmetry is broken, this theo- 
rem seems to indicate that T symmetry is also broken at 
the same level. These symmetry (and asymmetry) prop- 
erties have played an important role in developing the 
standard model of particle interactions, which describes 
electromagnetic, weak, and strong interactions. 

The discovery of CP- violation led to a fascinating ques- 



tion. Is it possible that particles, atoms and molecules do 
have permanent electric dipole moments after all? This 
would require interactions that violate both P and T, 
but we know that the weak and strong interactions to- 
gether can do that. As it happens, the standard model, 
with its standard P and T violation, predicts exceedingly 
small EDM values. This is the result of a fortuitous can- 
cellation, which comes about from the simplicity of the 
standard model. The pressing issue today is to discover 
what lies beyond the standard model. In order to un- 
derstand more clearly the origin of mass and in order 
to accommodate a quantum theory of gravity it seems 
very likely that there are more particles than the stan- 
dard collection. Such an increase in complexity immedi- 
ately leads to the prediction that particles, atoms and 
molecules have EDMs much larger than the standard 
model values. They are still very small, but they are 
no longer too small to measure. Therefore, the search 
for a permanent electric dipole moment of an atom or 
molecule is really the search for particle physics beyond 
the standard model. Molecules are beginning to play a 
central part in this search. 

The EDM de.p.nS of an electron, proton or neutron, 
is neccessarily aligned along the spin direction a of the 
particle. In essence, an EDM measurement in an atom 
or molecule involves polarizing the system with an ap- 
plied external electric field and searching for the interac- 
tion rjdx(J ■ E between the electronic or nuclear EDM and 
the polarised atom/molecule. Schiff's theorem [18] states 
that 7/ = if the atom/molecule is made of point parti- 
cles bound by electrostatic forces. In other words, the 
electronic or nuclear EDM does not see the applied field 
because it is shielded out by the other charged particles. 
This theorem is important for its loopholes: nuclei are 
not point particles and the electric dipole interaction is 
not screened when the electrons are relativistic. Conse- 
quently, rj is not zero if the atom/molecule is well chosen 
[TOl [20] . For example the best measurement of the pro- 
ton EDM comes from a measurement on TIF molecules 
[21] , where the large size of the Tl nucleus ends up giving 
7] ~ 1 for the nuclear spin EDM interaction. The upper 
limit on the neutron EDM is known both directly, from 
measurements on free neutrons [22j, and indirectly, from 
nuclear spin measurements on Hg atoms [23]. 

For electron EDM measurements, as opposed to neu- 
tron or proton measurements, rj can be much larger than 
1 if the electron moves relativistically within the atom or 
molecule. For example, the Tl atom is very sensitive to 
the electron EDM c?e, with 77 — —585. Currently, the best 
limit on the electron EDM is derived from such an exper- 
iment on Tl atoms [23]. The effective interaction -qdeS-E 
is linear in E because the polarisation of the atom is pro- 
portional to the applied electric field. This polarisation 
is small for atoms in laboratory strength fields because 
it derives from the mixing of higher electronic states in- 
duced by the field and these are typically 10^^ Hz higher 
in energy. By contrast, a heavy polar molecule is po- 
larised by mixing rotational states, which are typically 
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only 10 GHz away, giving five orders of magnitude more 
polarisability and a correspondingly larger rj. The po- 
larisation due to rotational mixing stops increasing once 
the molecule is largely aligned with the field and then 
rjE approaches a saturated value. For the YbF molecule, 
this value is an enormous 26 GV/cm [25] , 

A group at Imperial College, including the authors of 
this article, are currently in the process of measuring c?e 
using the YbF molecule [55] . In comparison with the Tl 
atom, this molecule gives roughly 500 times more EDM 
interaction energy, whereas the interaction with magnetic 
fields is essentially the same. Since stray magnetic fields 
constitute the primary source of systematic errors this is 
a significant advantage. However, Tl beams have much 
better statistical noise because they are much more in- 
tense than YbF beams. At present the gain in sensitivity 
is roughly offset by the loss in signal, and the YbF experi- 
ment is taking data at a level of precision similar to the Tl 
experiment [26j . In the new era of high-precision molecu- 
lar beam measurements that we are discussing here, the 
need for brighter, colder, slower sources is a recurring 
theme that we address again in the next section. 

Other molecular approaches to EDM measurement are 
also being pursued. The group of DeMille is aiming to 
measure the electron EDM in a metastable fi-doublet of 
PbO [57] using vapour in a cell. The group of Cornell is 
investigating the possibility of using trapped molecular 
ions [28 , among which HfF+ appears to be a promising 
candidate. There is also a proposal to make a dense 
sample of some radical such as YbF and to measure the 
magnetisation induced by aligning the electron EDMs in 
an applied electric field [29] . 

We turn now to parity violation without T violation. 
This is very well understood in the standard model of 
particle physics as a normal feature of weak interactions. 
It is also well established in atomic physics through the 
measurements on Cs atoms |30j . However, it remains 
a fascinating topic in the context of molecular physics, 
partly because it has not been observed in molecules but 
mainly because chirality plays such an important role 
in chemistry. The weak interaction is predicted to al- 
ter the energy spectrum between enantiomers of chiral 
molecules. Indeed, it is still a subject of debate whether 
the parity violation of weak interactions has played any 
role in establishing the chirality of the biochemistry in 
living organisms |3H 132] . The most intensively studied 
species are the methyl halides, CH-XYZ, where X, Y, Z 
are three different halogens [33] . The largest effect is pre- 
dicted to be a 50.8 mHz shift in the C-F stretching mode 
between left and right handed versions of CHFBrl. By 
contrast, the best experimental results [53] reach a preci- 
sion of 50Hz. It appears that cold trapped molecules will 
be necessary to measure weak interactions in these sys- 
tems. Much larger enantiomer shifts, in the range of sev- 
eral Hz, have recently been predicted for some rhenium 
and osmium complexes [3S] and these may be observable 
using supersonic beams. 

There is also nuclear physics interest in parity violation 



because it plays a role in nuclear structure. A particu- 
larly interesting possibility is that weak interactions in 
nuclei can induce an anapole moment, a P-odd multipole 
that produces no external field and corresponds in low- 
est order to a toroidal flow of current within the nucleus. 
An experiment is underway at Yale University to mea- 
sure the anapole moment of the ^"^^Ba nucleus using BaF 
molecules [35]. 

The last symmetry we mention here is Lorentz invari- 
ance, which has been a central plank of 20th century 
physics and has so far shown no indication of being vio- 
lated in nature. Even so, it is possible that new physics, 
associated with quantum gravity at the Planck energy 
scale, could lead to very small violations of Lorentz in- 
variance in the laboratory |37| . For example, there might 
be a change in the energy of an atom depending on the 
orientation of its spin relative to some preferred direc- 
tion, such as its velocity in the rest frame of the universe. 
Very sensitive experiments of this sort have already been 
performed using a variety of atomic clocks [38]. More 
recently, it has been pointed out that diatomic molecules 
provide a new way to investigate Lorentz invariance by 
orienting the internuclear axis of the molecule relative to 
the proposed preferred direction. The symmetry viola- 
tion could then be read out as a shift in the energy, bond 
length, vibration frequency, or rotation frequency |39j . 
Sensitivities for H2 and HD and their cations have been 
calculated in [33] • The authors conclude that an exper- 
iment to measure the ground state rovibrational transi- 
tions can improve limits on some elements of the Lorentz 
tensor c^^ [37] |40] by an order of magnitude. It is very 
likely that other more polar molecules could be conve- 
nient to use and that this area will develop as techniques 
for preparing, cooling and trapping molecules progress 
further. 

In the rest of this chapter, we discuss the production 
and detection of intense molecular beams, particularly of 
heavy polar radicals. We also describe methods of using 
pulsed molecular beams to map fields in the beamline 
and to detect interactions through the quantum coher- 
ences. Finally, we describe advances in slowing heavy po- 
lar molecules with a view to trapping them. These rather 
practical issues are the key points to be addressed if 
molecules are to fulfill the fantastic promise that we have 
just outlined for elucidating exotic fundamental physics. 



III. BEAMS OF COLD POLAR RADICALS 

All molecular beams must begin with a source. A 
good source will provide a large number of molecules 
in the quantum state of interest for the experiment. In 
most cases, the molecules need to be prepared in a single 
low-lying rotational state, and so the temperature of the 
molecules should ideally be smaller than the rotational 
energy spacing, typically IK or less. Pulsed sources of 
cold molecules, very narrowly distributed in both posi- 
tion and velocity, offer many advantages. They can be 
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prepared with very high intensity without imposing an 
excessive gas load on the vacuum system, they can be 
used to map the electromagnetic field along the beamline 
with high resolution and precision, they allow quantum 
coherences to be prepared and manipulated with great 
accuracy, and they can be slowed down in a Stark decel- 
erator to increase the coherence times. In this section, 
we concentrate on the formation and detection of cold, 
pulsed beams of the radical molecules that are typically 
required for measuring the P and T violating interactions 
discussed above. 

Supersonic expansion is a very common technique for 
producing cold molecular beams [H] . A high pressure 
gas expands through a nozzle into a vacuum chamber, 
acquiring a high centre-of-mass velocity but very narrow 
velocity distribution. The translational degrees of free- 
dom are cooled, as are the internal degrees of freedom of 
the molecules. While the first beams were continuous, the 
method was later extended by using pulsed valves with 
short opening times [13J |31] • The individual pulses could 
then be made very intense, without the gas load becom- 
ing excessive. Often, the molecules of interest have low 
vapour pressure, and need to be formed by laser ablation 
or electric discharge techniques. This is usually done im- 
mediately outside the nozzle of the pulsed valve, or inside 
an extended nozzle, where the density of carrier gas is 
high enough to entrain a useful fraction of the molecules 
produced |1S]. In our laboratory, we have used laser ab- 
lation to produce cold beams of YbF |46] , CaF and LiH 
molecules [17|- In all cases, we detect our pulsed beams 
using time-resolved, Doppler-free, laser induced fluores- 
cence (LIF). This detection method is very well suited 
for precision measurements where high sensitivity, high 
frequency resolution, and good beam diagnostics are all 
required. 



A. Apparatus 

The typical experimental setup is illustrated in Fig. IT] 
A solenoid valve emits short pulses of the carrier gas [48] . 
which is usually Ar, Kr or Xe, into a vacuum chamber 
maintained at pressures below 10~* mbar. The central 
part of the gas pulse passes through a skimmer, of diame- 
ter 1-2 mm, placed 50-100 mm downstream of the source, 
and into the high vacuum region where the pressure is 
below lO"'^ mbar. A fast ionization gauge placed on the 
axis of the beamline can be used to ensure good align- 
ment. Using two such gauges, one immediately outside 
the nozzle and the other much further downstream, the 
initial width, the speed and the translational tempera- 
ture of the gas pulses can all be measured. The shortest 
gas pulses we obtain with this valve have a full width at 
half maximum (FWHM) of 81 /iS 135]. In an ideal super- 
sonic expansion from a reservoir at temperature Tq, the 
carrier gas approaches a terminal velocity 
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FIG. 1: Typical experimental setup for producing pulsed 
beams of cold polar radicals, and detecting them by laser 
induced fluorescence. Not to scale. 



where 7 is the specific heat ratio, 5/3 for an ideal 
monatomic gas, and m is the mass of an atom of the 
carrier gas. Our measurements using a range of carrier 
gases and temperatures show that the true speed of the 
carrier gas is 5-15% faster than vt- 

A suitable target, placed outside the nozzle of the 
valve, is ablated using light from a Q-switched laser. To 
ensure minimal disruption of the gas pulse, the target is 
thin in the direction of the gas-jet, typically 2-5 mm. The 
resulting ablation plume can usually be observed by eye. 
It is forward peaked along the normal to the target, which 
is usually perpendicular to the direction of the beamline. 
Some of the atoms and molecules in the plume become 
entrained in the high density carrier gas, thereby cool- 
ing to temperatures and speeds that approach those of 
the carrier gas. The target may either contain all the 
precursors needed to form the molecules of interest, or 
some of the precursors can be added to the carrier gas. 
For example, we have produced beams of YbF cither by 
ablating pure Yb and mixing SFg into the carrier gas, 
or by ablating a pressed target containing a mixture of 
Yb and AIF3 powders. The former method is slightly 
simpler to realise, but both methods work equally well. 
We typically use ablation pulses of 5-10 ns duration, fo- 
cussed to a spot-size of 1-2 mm. Under these conditions, 
the optimal ablation energy is in the range 10-50 mJ. 

A photomultiplier, placed some distance (typically 10- 
150 cm) downstream, detects the molecules by means of 
cw-laser-induced fluorescence. The laser is directed at 
right angles to the molecular beam in order to minimise 
the Doppler shift and is tuned to a strong molecular 
transition. The transit time of the molecules through 
the laser beam is typically 5 — 10 ^s, giving this method 
high temporal resolution. The width of the fluorescence 
excitation spectrum is usually greater than the transit 
time limit, being typically 20-50 MHz. This comes from 
a combination of the residual Doppler width, due to the 
angular distribution of the molecules, and the natural 
lifetime of the excited state. Such high spectral resolu- 
tion is typically needed in order to resolve the hyperfine 



5 



structure and hence to permit coherent state manipu- 
lation and readout at rf frequencies. The fluorescence 
detection method is very sensitive. If the detection effi- 
ciency is e (assumed <C 1) and there are N molecules per 
shot passing through the detector area in a time interval 
w' , then the shot noise limited signal-to-noise ratio for a 
single shot will he s : n — eN/ \/bw' -\- eN, where b is the 
background rate of detected photons. In a good detec- 
tion setup, bw' ~ 1 and e is in the range 1-10%. Taking 
e = 0.02 as a typical value, we find that the signal-to- 
noise ratio per shot is 1 when = 81 molecules. 



B. Translational temperature and source size 

Figure |2] shows the time-of-ffight profiles of ground 
state YbF molecules recorded by two separate LIF detec- 
tors placed 340 mm and 1300 mm from the pulsed source. 
The background laser scatter is lower at the downstream 
detector. In this experiment, the molecules were pro- 
duced by ablating pure Yb just outside the nozzle of 
the solenoid valve which was pressurized to 4 bar with 
a mixture of Ar(98%) and SF6(2%). The probe laser was 
tuned to excite the F = 1 hyperfine component of the 
X'^^+{v" 0) - A'^Ily2{v' = 0) Q(0) transition. There 
are four contributions to the spread of arrival times - (i) 
the temporal spread of the molecules produced at the 
source, (ii) their spatial spread at the source, (iii) the 
forward velocity distribution in the pulse and (iv) the 
temporal resolution of the detector. The last of these is 
usually small enough to be neglected. 

The fiux of molecules with speeds in the interval 
from V to V + dv is usually taken to be f{v) dv — 
Av^ exp{—M{v — vo)"^ /^ksT) dv, where M is the mass, 
T is the translational temperature, vq is the central ve- 
locity and A is a normalizing constant. Consider those 
molecules born in the source at time tg with initial posi- 
tion s along the beam axis. In the LIF detector, placed a 
distance L away from the source, these molecules produce 
a time-dependent signal 



h{t,ts,s) 



-Aftig {to-s/vo-t+t.y^ 
2kBT {t-t,y 
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where to — L/vq. As we shall soon see, both the temporal 
width and spatial width of the source are small, ts <C to 
and s ^ L. Furthermore, at typical detection distances, 
the range of arrival times is much smaller than the mean 
arrival time, and it is valid to set t « everywhere 
except in the numerator of the exponent. With these 
approximations, Eq. ([2| simplifies to 
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FIG. 2: Time of flight profiles of ground state YbF recorded 
by two laser induced fluorescence detectors, situated 340 mm 
and 1300 mm from the source. The lines are Gaussian fits to 
the two datasets. 



The signal at the detector is obtained by integrating 
over the temporal and spatial distributions present in the 
source. We do not know these distributions, but we can 
hope to obtain a measure of their characteristic widths. 
In this spirit, we assume that the source emits the nor- 
malised distribution 
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Then, the detector records the signal 
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where the pulse width w 
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broadening due to the distribution of positions and times 
where molecules are first formed. 

As shown by the Gaussian fits in Fig.|2] the recorded 
profiles are well described by the model except in the high 
velocity tails. The long tail, indicating a hotter compo- 
nent in the beam, is a common feature of our source 
when optimized for maximum signal and minimum shot- 
to-shot fluctuation; it can be removed by re-optimizing 
for low temperature. By using two well-separated detec- 
tors, the translational temperature can be obtained. It 
is 



\~A^^ ( 41n2(f-V^i^p)2\ 
h{t,ts,p) K A-^ey.^\^ j, (3) 

1 /2 

where p — —s/vn, and w = [81n2 ksTtl/ [Mvl)) is 
the temporal width (FWHM) of the pulse due to the 
thermal spread of forward velocities. 



81n2fcB tl~t\ ' ^ ^ 

where w'-y, w'^, i\ and ti are the widths (FWHM) and 
central arrival times obtained from Gaussian flts to the 
downstream and upstream data. For the data shown, the 
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speed is vq = 586 m/s and the temperature is T = 4.8 K. 
The use of two detectors allows an unambiguous deter- 
mination of the temperature, whereas a single detector 
alone sets only an upper limit. However, if the detector 
is far from the source, that upper limit can be very close 
to the true temperature. For this data, the downstream 
detector alone provides an upper limit that is just 1.3% 
higher than the measured 4.8 K. 

Having measured the temperature, we can also extract 
from the data upper limits on the initial temporal and 
spatial spreads at the source. For the data shown, these 
are At^ < 14.5 /iS and voAp < 8.5 mm. In separate ex- 
periments, using a dual ablation technique, we have mea- 
sured At^ « 5 /zs H?] . 



C. Molecular flux 

We consider next how to determine the absolute flux 
of molecules from the LIF signal. The detector counts 
the flux of photons, and we wish to convert this to a 
flux of molecules by knowing the mean number of flu- 
orescent photons emitted from each molecule. To flnd 
this, we model the molecule as a 3-level system and use 
rate equations. Every molecule starts out in level 1 and 
passes through the laser, which excites the resonance be- 
tween levels 1 and 2. The excitation rate is R, as is the 
rate of stimulated emission. It is directly proportional 
to the laser intensity, /, and depends on the detuning 
S — lol —uJi2 of the laser angular frequency, lol, from the 
molecular resonance frequency, ^12. Taking the optical 
Bloch equations in the limit where the coherences have 
reached a steady state, one finds that 



(i + 452/r2)// ^ ' 

where F is the spontaneous decay rate of level 2, = 
tQch^ir 12)^ / is the saturation intensity, and D is the 
matrix element of the dipole operator connecting levels 1 
and 2. Level 1 is stable, while level 2 decays with rate rV 
to level 1, and rate (1 — r)r to level 3 which represents all 
the other states in the molecule. The rates for excitation 
or decay out of level 3 are negligible. 

We solve the rate equations to find the number of 
molecules in level 2 as a function of time, N2{t). Inte- 
grating TN2{t) over the laser-molecule interaction time, 
r, we find the number of fluorescent photons emitted per 
molecule to be 

^ Rjf^ — R— \ R— J 

where 

R± ^ R + T/2±y/R^ + rRT + T^I'i. (9) 




FIG. 3: Number of fluorescent photons per molecule, Up, 
when the branching ratio for returning to the initial state is 
r = 0.5. In the main graph, Up is plotted versus the intensity, 
I / Is, for T T = 10, 20, 50, 100. The inset shows for each case 
how Up varies with the laser detuning, 5/V, when I — Is/3. 

In the limit where i?+r 3> 1 and i?_T 3> 1, ?^p acquires 
its asymptotic value np. max = 1/(1^'')! this is simply the 
sum of the obvious geometric series, rip^max = X]??=o^^- 
For good detection efficiency, we would like to ensure that 
Up reaches its maximum possible value, and so we should 
examine how easily this limit is reached. If, as is usual, 
level 2 is an electronically excited state with an allowed 
electric dipole transition, F will typically exceed 10^ s~^. 
The interaction time, r, is usually greater than 1 fis, and 
so we are in the limit Ft 3> 1. If we suppose that the 
laser excitation is weak, in the sense that i? <C F (or, 
equivalently, / <C /s), we get 

1 - e-«(i-'-)^ 

np= (FT>l,i?<F). (10) 

1 — r 

The limit Up = 1/(1 — r) is reached once i?(l — r)r ^ 1, at 
which point the fluorescence signal saturates. Note that 
it is usual for Fr to be greater than 100, and that for 
molecules, it is rare to have r close to 1. It follows that 
the above 'saturation' condition can easily be met, even 
when I <^ Is - the fluorescence signal saturates for laser 
intensities well below Is because the interaction time is 
very long compared to the time required for the molecule 
to reach level 3, the dark state. 

Figure [3] shows the value of Tip clS function of I /Is for 
four different values of Ft with r — 0.5. The maximum 
value of Up is 2, and in the case of Ft — 100, Up is 
within 1% of this limit at / = /s/4. At lower values of 
Ft, the asymptotic value is lowered and more intensity 
is required to reach it. For example, when Ft = 10, Up 
has an asymptotic value of 1.84 and is 78% of this value 
when I — Is- The saturation process described here also 
leads to a type of power broadening of the spectral line 
when the laser frequency is scanned. An increase in laser 
intensity increases the fluorescence in the wings of the 
resonance more than at the centre of the resonance, and 
so the line is broadened. This is seen in the inset to Fig. [3] 
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where the lineshape is plotted for the various values of 
Ft, in the case where r = 0.5 and / = /s/3. 

If the decay rate of the excited state happens to be 
small, it is possible that the system will be in the opposite 
limit, Ft <C 1. In that case, even if ^ F, the mean 
number of photons scattered by each molecule will be 
much smaller than 1, and independent of r. We cannot 
then use Eq. ([8| because the damping has no time to 
act. However, Eq. (20l, which has no damping, can be 
integrated to give 




sm{TT^I/{2Is)) 



(Ft< l,i?>F). 

(11) 

The fluorescence saturates to the value Ft/2 at high laser 
intensities. 

Finally, it is interesting to consider the case where the 
interaction time is long {Rt 1, Ft ^ 1) and r is 
very close to 1 so that each molecule has the possibility 
of scattering a large number of photons. Then, Eq. ([s]) 
reduces to the familiar atomic physics result, 



Ft 



I/Is 



RTt 



2R + T 2 1 + ///, + 4(57F2 



(12) 



The on-resonance fluorescence now saturates to Ft/2 
when the condition I ^ Is is met. 

Returning to the determination of the molecular flux, 
we write the total number of photons detected per shot 
as 



Eq. Q and then using Eqs. ^ and ^ (the probe laser, 
of wavelength A, is assumed to be on resonance). The 
calculation requires some knowledge of r, F and Ig, but 
provided the detection is in the saturated regime, rough 
estimates will suffice since the result becomes rather in- 
sensitive to F and Ig, and also to r if r is small, as is 
often the case. 

For the detection setup shown in Fig.[l] the total de- 
tection efficiency is 



e = (f^i/4^) J2 9.(1 + n{Km{\)\{K) . (15) 



Here, qi is the fraction of fluorescent photons in the emis- 
sion line whose wavelength is A^, TZ, Ti and x &re the 
wavelength-dependent mirror reflectivity, lens transmis- 
sion, and photodetector quantum efficiency and f2i is the 
solid angle subtended by the light-gathering lens. If win- 
dows and filters are present in the setup, their transmis- 
sions need to be included too. 

By measuring the value of p and performing the above 
calculations, the flux of molecules in the detected quan- 
tum state can be determined. With careful measure- 
ments, an uncertainty below 50% should be possible. For 
the cold YbF molecules produced in our laboratory, the 
flux is measured to be 1.4 x 10^ ground state molecules 
per steradian per shot, when the carrier gas is argon |46j . 



D. Rotational temperature 



L2 



N{x, z)e(x, z)np(x, z) dx dz, (13) 



where N(x,z), e{x,z) and np{x,z) are the number of 
molecules per steradian in the quantum state of interest, 
the detection efficiency and the number of fluorescent 
photons per molecule, all in the interval dx dz around 
the point (x, z) in the plane of detection. To determine 
the molecular flux as accurately as possible, the detec- 
tion area should be defined by placing a small aperture 
on the molecular beam axis, slightly upstream of the de- 
tector, to block all but a small portion of the molecular 
beam. Similarly, the probe laser should be coUimated 
and prepared with a top-hat intensity distribution in the 
X direction. Then N and e will be constant over the de- 
tection area. The value of Up will be independent of x, 
but because of the Doppler shift, it remains a sensitive 
function of z, particularly if the degree of saturation is 
not high. Under these conditions, Eq. ( 13 1 reduces to 



HNe 
L 



np{e) d0, 



(14) 



where 9 — z/L and H is the height (along x) of the 
detection area. The integral is straightforward to calcu- 
late by substituting the Doppler shift 6 = 2'kvqOIX into 



The rotational temperature of the molecules can be de- 
termined by scanning the laser frequency and recording 
the rotational spectrum. When the rotational tempera- 
ture is Tr, the intensity of a rotational line in the spec- 
trum is proportional to N{J)e{J)np{J), where N{J) — 
{2J+ 1) exp[—BJ{J+ l)/kBTr] is the relative number of 
molecules in the rotational state J, B is the rotational 
constant, and e(J), np{J) express the J-dependence of 
the detection efficiency and the number of scattered pho- 
tons per molecule. Usually, the efficiency e is almost 
independent of J. In many cases this is also true of Up, 
making it easy to extract the rotational temperature from 
the relative line intensities once a few rotational lines 
have been measured. A more accurate temperature de- 
termination must take into account the variation of the 
matrix elements with J. For example, in — ^ S tran- 
sitions, the Mj-averaged value of in the R-lines is 
proportional to {J+ l)/(2J-t- 1), which decreases from 1 
to 1/2 as J goes from to oo. For the P-lines, by con- 
trast, cx J/{2J + 1), which increases from to 1/2. 
Since this variation of the matrix elements generally af- 
fects both the excitation rate R and the branching ratio 
r, it influences the value of np(J) in both the saturated 
and unsaturated regimes. In the saturated regime, where 
rip ~ 1/(1 — r), the variation of Hp with J is strongest if 
r happens to be close to 1. 
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In our lab, we have measured rotational temperatures 
of YbF, CaF and LiH beams. For the first two, the rota- 
tional temperatures are usually very close to the trans- 
lational temperatures (typically in the range 1-5 K) 
For LiH we measure rotational temperatures consider- 
ably higher than the translational temperature [47] , 



E. Source noise 

The flux of molecules obtained from these sources is 
subject to shot-to-shot fluctuations, as well as having a 
slow drift (mostly downward) . The time scale for the slow 
drift is typically 10* — 10^ shots on a given spot of the 
target. We usually attach the target to the rim of a large 
disk, typically 20 cm in diameter, which we rotate incre- 
mentally as each target spot becomes exhausted. In this 
way, the lifetime of a target is of order 10^ shots. The 
ability to run the source continuously for long periods 
of time is very important for precision measurements. A 
second essential requirement is that the shot-to-shot fluc- 
tuations be small, since they can contribute directly to 
the noise in the experiment. In measuring the electron 
electric dipole moment using cold YbF molecules, our 
detector records approximately 3000 photons per shot, 
with a corresponding \/]V photon shot-noise limit of 2%. 
Source fluctuations should ideally be kept below that 
level, but this proves rather difficult to achieve; when 
optimized, the short-time-scale fluctuations of our source 
are typically 2-3%. 



IV. COHERENT MANIPULATION OF 
INTERNAL STATES 

A. Stark and Zeeman shifts of the hyperfine states 



In section |III C| we discussed the excitation of higher 
electronic states using laser light to drive optical dipole 
transitions. The interaction was strongly damped by 
spontaneous emission from the upper level; indeed, the 
molecules were detected by means of the scattered pho- 
tons. In this section, we consider some ways to manip- 
ulate the hyperfine sublevels within the ground state of 
the molecule. In contrast to optical transitions, the co- 
herences between these ground-state levels are not radia- 
tively damped because the spontaneous transition rates 
are very low for transition frequencies in the sub-GHz 
range. The control of these coherences provides a basis 
for exceedingly high precision measurements of electric 
and magnetic fields and for measurements such as that 
of the electron EDM. 

In order to be concrete in our discussion, let us take 
the simple case illustrated in Fig. |4] of two hyperfine lev- 
els, F — 1 and F = 0, such as one finds in a diatomic 
X^S(A^ = 0) molecule with nuclear spins of and 1/2. 
Here the total angular momentum F is the sum of the 
electronic angular momentum J = 1/2 and the nuclear 



F=1 



(1,0) 
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FIG. 4: Hyperfine levels. Left: field-free levels F = 0,1. 
Centre: electric field induced shift of levels and tensor Stark 
splitting of the triplet F — 1. Right: Zeeman splitting of the 
doublet F = 1, Mf = ±1. 



spin / = 1/2. Two examples of such molecules, which 
have been studied in our laboratory, are ^^''YbF and CaF. 
In the absence of any external fields, the two hyperfine 
levels are separated by the hyperfine splitting A, as il- 
lustrated on the left of Fig. [4] and the three magnetic 
sublevels of = 1 are degenerate. {A ~ 2tt x 123 MHz 
in CaF and 2tt x 170 MHz in YbF). For a fuller discus- 
sion of hyperfine structure in such molecules, see |49j . 
When the molecule is subjected to an electric field E^, 
the main effect is the rigid rotor Stark shift due to the 
electric dipole moment fi^ along the internuclear axis of 
the molecule, as we discuss more fully in Sec|V] This shift 
is large in comparison with the hyperfine interaction; for 
example the = state of YbF shifts downwards by 
20 GHz at 20kV/cm. To a good approximation all four 
hyperfine levels shift together, but there are some small 
differential shifts as well, due mainly to the tensor part of 
the hyperfine interaction between J and /, the electronic 
and nuclear angular momenta. This effect is analysed in 
detail in Ref. [13]. Relative to the F = level, the state 
{F,Mp) = (1,0) shifts up by Aq whereas states (1,±1) 




5 10 15 
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FIG. 5: Stark shift of the ground state hyperfine intervals in 
YbF. 
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shift up by Ai, as illustrated in the centre of Fig. [4] 
These shifts are plotted in Fig. [5] versus electric field for 
the particular case of YbF. The lines are a calculation 
following the theory detailed in , which has been con- 
firmed by experiment. One sees that these shifts of the 
hyperfine frequencies within the = manifold are typ- 
ically a thousand times smaller than the over-all shift of 
the manifold itself. The general behaviour is a quadratic 
shift at low electric fields, where the characteristic dipole 
interaction — /Xe-E is small compared with the rotational 
constant B, evolving to a more nearly linear behaviour 
at higher fields and ultimately saturating at the highest 
fields. 

It is worth noting that the total Hamiltonian (includ- 
ing the coupling to the external electric field) is invariant 
under time reversal. Since the two states {(1, -|-1) in field 
E} and {(1, —1) in field E} are time-reverses of each other 
and since the Hamiltonian is invariant under time rever- 
sal, it follows that the energy levels of (1, ±1) are exactly 
degenerate at all electric fields. Thus, even though the 
molecule has an electric dipole moment along its inter- 
nuclear axis, and has an induced electric dipole moment 
along the applied field direction z, it does not have an 
electric dipole moment proportional to F^. This is a di- 
rect consequence of time reversal symmetry. By contrast, 
if the electron were to have a permanent EDM dg along 
its spin, this would lift the degeneracy between the two 
levels - a direct consequence of the fact that such an EDM 
violates time reversal symmetry. This spin-dependent 
Stark shift is what our group measures in searching for 
an electron EDM using YbF molecules [23]. 

A more mundane way to split the (1, ±1) levels, as il- 
lustrated on the right of Fig. |4] is with a static magnetic 
field B through the Zeeman interaction g^/i^E-B, where 
/is is the Bohr magneton, and the g- factor expresses the 
ratio of magnetic moment to total angular momentum. 
Provided this interaction is small compared with the ten- 
sor Stark splitting Ai — Aq, the molecule only responds 
to the component parallel (or antiparallel) to the elec- 
tric field. This splits the levels by iA^ = zLgp^sBz/fi. 
The additional contribution to this splitting resulting 
from a perpendicular field component B± is of order 
^zifJ-sB^/iAi — Aq)]^ and is therefore negligible. 
No violation of time reversal symmetry is implied here 
because the state {(1,-1) in field B} is not the time- 
reverse of {(1, -f 1) in field B}: time reversal also reverses 
the sign of B. There is no shift (to first order) of the states 
(F, Mf) = (0, 0) and (1,0) since they have < F, >= 0. 

To summarise, the mean shift of the two levels (1,±1) 
relative to level (0, 0) is a measure of the electric field 
strength, whilst the splitting between the two measures 
the parallel component of magnetic field (and a possi- 
ble very small contribution from the electron EDM) . Al- 
though we have chosen to illustrate this with a simple 
hyperfine system having F — and F = 1, the same 
behaviour applies more generally, namely that any two 
levels F, zLMp are equally shifted by the Stark interac- 
tion with electric field Ez, whereas they are split apart 



by a magnetic field along z (and by an EDM). In the 
coherent manipulations we discuss below, the (1,0) level 
plays no significant role and we therefore do not consider 
it further. The remaining three levels can now be conve- 
niently abbreviated to (0), (-1-1) and (—1), with energies 
0, A + Ai + Az and^-|-Ai— A^ respectively. In the next 
section, we will find it convenient to redefine the zero of 
energy. 



B. Two-pulse interferometry of a 3-level system 

Let us write the amplitudes of these three states as a 
column vector 



ao{t) 
az{t) = I a+i{t) 
a-i{t) 



(16) 



where the subscript z indicates that the quantisation axis 
is taken to be along the electric field. The free evolution 
of these amplitudes from time ti to time -I- r in the 
presence of the static electric and magnetic fields is given 
by the propagator 



no(ti,T) = 








(17) 



such that az{ti + r) = no(ti, T)az(ii). Here the Stark- 
shifted hyperfine interval A -|- Ai has been replaced by 
the symbol CI and we have moved the zero of energy to 
the centre of that interval in order to simplify the algebra 
that follows. 

Let us suppose that the molecules are prepared in state 
(0), then subjected to an rf magnetic field cos{ujt+(j)) 
along X in order to drive transitions to states (-1-1) and 
(—1). This field excites the coherent superposition state 
(c) = ^[(+1) + (^1)] and docs not couple at all to the 
orthogonal superposition (u) = -i=[(-|-l) — (—1)]. (The 
converse is true for an rf field along y). This suggests a 
new basis with quantisation along the x-axis, in which 
the new state amplitudes are given by 




UUz 



(18) 




Note that the transformation U is its own inverse: 
U = U^^. In the x-basis, only the states (0) and (c) 
are coupled by the rf field. With this reduction of the 
problem to a two-level problem, we can write down how 
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the amplitudes evolve in the a;-basis under the influence 
of the rf magnetic field applied from time ti to time ti+r. 
Following the standard derivation leading to Eq. V.7 of 
Ramsey's book [50], we find that 



(19) 



e"'?^ 



where 



Z — jcosHsm( — ) + cos( — ) 
W — isin0sin( — ) 



cos 6 = 



le = - 



a 
2b 



fJ-x^x |c) 



Although the phase of the field (f> is not important when a 
single pulse is applied, we keep it in this formula because 
it becomes relevant when we consider double pulses. As 
the state (c) is excited, the effect of any static magnetic 
field Bz will be to rotate it into the third state (u) at the 
Larmor frequency A^. This has been ignored in deriv- 
ing Eq. (19), under the assumption that the rf excitation 



will be performed quickly in comparison with the Lar- 
mor precession. Taken together, Eqs. (17), (18) and (19) 



provide us with all the tools we need to investigate the 
evolution of this three-level system under any sequence 
of short rf pulses. 

When molecules in state (0) are subjected to a single 
rf pulse, the probabihty of excitation to state (c) is given 
by Eq.ra as 



^1 pul. 



.(0- 

46^ 



w)2-K4fe2). (20) 



This is the usual magnetic resonance lineshape for 
transitions in a 2-level system without damping. At res- 
onance the population oscillates sinusoidally between the 
two states (this is known as Rabi oscillation). A 'vr-pulse' 
is an on-resonance pulse with 26t = tt, which transfers 
all the population from state (0) to state (c). In section 
II V CI we will discuss how this can be used in a molecular 
beam to map out the fields along the beamline. An on- 
resonance '7r/2-pulse' (26t = 7r/2) drives the transition 
only half way, creating an equal superposition of states 
(0) and (c) with a definite relative phase. The density 
matrix element describing this coherence at the end of the 
pulse (at time ti+r) is {aoac*)x — \i exp{i[ti;(ti-|-T)-|-</)]}, 
the phase of which (apart from the fixed factor of i) is 
just the final phase of the rf field. 



In conventional Ramsey spectroscopy of a 2-level sys- 
tem [50], two short 7r/2 pulses are applied in succession. 
If the second pulse comes immediately after the first, the 
transition is completed and all the population is excited. 
If instead there is a delay time T between the two pulses, 
which is long compared to the pulse duration r, the tran- 
sition probability becomes sensitive to small differences 
between the rf frequency and the molecular transition 
frequency. The internal coherence evolving at the tran- 
sition frequency accumulates a phase between pulses of 
whereas the rf field evolves a phase ujT. When the 
difference between these two reaches tt, the second pulse 
reverses the effect of the first, returning all the popula- 
tion to the initial state. More generally, the probability 
that a molecule will end up in the excited state is 

PRa,nsey{Q ^ 1) = ^{1 + COs([r! - Uj]T - <50)} , (21) 



where we have allowed for the useful possibility of giving 
the second pulse a phase 4> + 6(j) when the first has phase 
(j). These oscillations of the population, resulting from 
a beat between the coherence and the driving field, are 
known as Ramsey fringes. They are important because 
a long waiting time T allows a small difference f2 — w to 
be measured, giving rise to very precise spectroscopy. 

We have generalised this idea to our system of three 
hyperfine levels. The result, derived from Eqs. (17), (18) 
and (19) is 



P2xpulse{0 ^ 1) = (22) 

1 - cos2(5t)i cos2(6t)2 - cos2(A^r) sin2(6T)i sm^{bT)2 
+ ^ cos([r2 - uj]T - 6<l)) cos(A^T) sin(2&T)i sin(26T)2 . 

Here we allow for the possibility that the values of br for 
the two pulses, (6t)i and {bT)2, are not equal. In the 
'Ramsey' case of two 7r/2-pulses, this simplifies to 



i (3 - cos^iA^T) + 2 cos{[n - lu]T - 



(23) 



)cos(A^r)) 



If we set A2 to zero, Eq.(p3| reduces to the standard two- 



level Ramsey result of Eq.(21l because the third state 
(u) plays no role when the two levels (+1) and (—1) are 
degenerate. When A^ 7^ it can be useful to pick out the 
Ramsey interference by switching S(t> between and tt and 
taking the difference, which is cos([ri — uj]T) cos{ A ^T). 
As the frequency of the oscillator is swept, the amplitude 
of these fringes provides information about the Zeeman 
shift A^ , while the phase of the fringe pattern reveals the 
precise value of the splitting f2 with a precision controlled 
by the choice of T. In section |IVD| we give an example 
of how this can be used to look for small changes in a 
large electric field. 

With three levels, it becomes possible to do interfer- 
ometry using 7r-pulses as well. Taking 2bt = tt, Eq.(22l 
becomes 



2X7T-pulse 



(0-^ 1) =sin2(A3r) 



(24) 
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FIG. 6: Schematic diagram of the beam machine. Bunches 
of molecules issue from the source (in the y-direction) and 
are skimmed before being optically pumped into a single hy- 
perfine state. The molecules enter the magnetically shielded 
interaction region and fly through a high-voltage capacitor 
where electric and magnetic fields can be applied along z. 
This doubles as an rf transmission-line where the rf magnetic 
field is along x. Finally, the molecules are detected by laser 
induced fluorescence. 



This has a simple interpretation. Molecules are excited 
by the first pulse to state (c). This state subsequently 
evolves into cos(A2T)(c) -I- isin(A2T)(u) because of the 
splitting between states (+1) and (—1). The second pulse 
then drives the state-(c) part of the population back to 
(0), leaving those in state (u) alone. As the magnetic 
field is scanned, this produces fringes in the final state-(O) 
population, with a spacing that is inversely proportional 
to T . These interference fringes can be used for sensitive 
magnetometry and in searching for an electron EDM. 



C. Experiments with single pulses 

Precision measurements usually require careful con- 
trol and monitoring of stray and applied fields, both 
electric and magnetic, throughout the interaction region 
of the apparatus. The small spatial and temporal ex- 
tent of molecular beam pulses make it possible to do so 
with high spatial resolution [5T. Fig. |6] shows an ap- 
paratus to demonstrate this using YbF molecules and 
pulsed rf fields. A pump laser and an interaction re- 
gion have been added in between the source and detec- 
tor already shown in Fig.[l] Although the beam is cold, 
both hyperfine levels are occupied because the splitting 
A = 170 MHz is very much less than kT . The pump laser 
excites the A'^Ili/2 — X'^T,'^ Q(0) transition at 552 nm, 
for which the Doppler width is (almost) eliminated by 
pointing the laser beam perpendicular to the molecu- 
lar beam. This makes the excitation spectrum narrow 
enough (~ 20 MHz) to excite just the F = 1 population 
so that it becomes selectively depleted. The remaining 
= molecules are virtually all in the F = state, 
which serves as the initial state for subsequent manipu- 
lation by rf pulses. 

The interaction region starts 450 mm from the skim- 
mer and is 790 mm long. It is magnetically shielded to 
reduce the ambient field, whilst current-carrying wires 
inside the shields can generate a magnetic field if re- 
quired. Within this region there is a pair of electric field 
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FIG. 7: Excitation spectra for the ground state F = ^ 
F — 1 transition in YbF using 7r-pulses of three different 
pulse lengths (upper, t — 18 /xs; middle, t — 36 fia; lower, 
T — 72 /is). A static electric field of 12.5kV/cm is applied. 
Dots: experimental data. Lines: plots of Eq.|20[ 



plates, 750 mm long and 70 mm wide, with a 12 mm 
spacing, constant to better than 200 /im over the full 
length. These are machined from cast aluminium, then 
gold coated to improve the uniformity of the surface po- 
tential (using a non-magnetic, nickel-free process). The 
whole assembly is non-magnetic. With a field across the 
gap of 15 kV/cm the leakage current is less than 1 nA. 

The same plate structure also serves as a 34 trans- 
mission line for the rf field, transporting it as a 170 MHz 
TEM wave travelling parallel or antiparallel to the beam 
direction. This is described more fully in |51j . As a bunch 
of molecules travels along the interaction region, a hyper- 
fine transition can be induced at any desired position by 
pulsing the rf field on for a short time. This repopulates 
the F — I state and therefore produces an increase in 
the fluorescence signal at the detector, as shown by the 
rf frequency scans in Fig.[7l Being a TEM wave, the rf 
magnetic field between the plates is accurately perpen- 
dicular to the static electric field and we choose to define 
its direction as the x-axis. This field therefore drives the 
transition (0) — (c) discussed in the last section. 

Fig. [T] shows three excitation spectra, obtained by 
applying 7r-pulses of three different durations (18, 36, 
and 72 /zs) to molecules near the middle of the inter- 
action region. Superimposed on the data points are 
solid lines corresponding to Eq. ( 20 1 , which does a good 



job of describing the lineshapes, including the positions 
and relative heights of the sidebands. When we fix the 
value of 26r at some value 9 (equal to tt in this case), 
these lineshapes can be re-writen as a universal function 
sin2(|Vl + a;2)/(l + x^), where x = - u;)t/6'. Thus, 
the width of the line is inversely related to the duration 
of the pulse, becoming wider as the pulse is made shorter, 
as one would expect from the usual Fourier relation be- 
tween pulse duration and spectral width. For a 7r-pulse, 
the full width at half maximum is (5wfwhm = 5.0/t. 
The transition probabilities in Fig. [7] peak at 
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FIG. 8; Upper curve: the hyperfine transition moves to lower 
frequency when the electric field is reduced. The dashed line 
marks the line-centre of Fig.[7|at higher electric field. Lower 
curve: The line splits when a magnetic field is applied. 



173.513 MHz, not at the field-free frequency of 
170.254MHz, because these spectra were measured with 
a potential difference of 15 kV across the 12 mm gap of 
the plates and are therefore Stark- shifted by A i, as illus- 
trated in Fig. |4] The magnitude of this shift is a measure 
of the electric field strength at the place occupied by the 
molecules when the rf pulse was applied. By varying the 
timing of the rf pulse, it is possible to map out the elec- 
tric field as a function of position along the beamline [51]. 
Note that the Doppler shift of the rf transition is only 
a few hundred Hz and is therefore insignificant at this 
level of accuracy. In Fig. [8] we show how the resonance 
frequency moves down (by 40.246 kHz) when the applied 
potential difference is reduced (by 203 V). The spectral 
resolution afforded by the 18 /is 7r-pulses is quite suffi- 
cient to see this shift clearly. Fig. |8] also shows how the 
line splits by 2Az (see Fig. |4| when a DC magnetic field 
of 0.8 iiT is applied. In order to resolve the 22 kHz split- 
ting fully, the linewidth was narrowed by increasing the 
pulse duration to 72/is. Over this time, a YbF molecule 
travelling at 590 m/s covers 42 mm, so it has been neces- 
sary to give up some spatial resolution in order to achieve 
the higher spectral resolution. The use of long rf pulses 
for high precision spectroscopy is not ideal because the 
static field average that is measured can be affected in a 
complicated way by variations in the strength and polar- 
isation of the rf field, either in space or in time. As was 
first pointed out by Ramsey, it can be more satisfactory 
to use a pair of rf pulses, and that is what we discuss 
next. 



D. Experiments with double pulses 



In Eq. ( 23 1 , we found an expression for the lineshape 
using two short 7r/2- pulses separated by a time T. In 
order to pick out the Ramsey interference, we suggested 
introducing a phase shift between the first and sec- 





FIG. 9: Ramsey interference signals using two 7r/2-pulses sep- 
arated by a time T — 900 fis. The signals are plotted versus 
the rf frequency ui and also as a function of magnetic field in 
uT. (a) Theoretical signal, scaled to have the same amplitude 
as the experimental result, (b) Measured Ramsey interfer- 
ence, in excellent agreement with the theory. 



ond pulse and taking the difference between = tt 
and (j) — 0. This is expected to give fringes of the 
form cos(A2r) cos([ri — w]T), which we have plotted in 
Fig. [9]^ a) versus rf frequency and magnetic field for the 
case when T = 900 /is. In Fig. [9|b) we display the re- 
sult of an experiment to test this formula, in which we 
scanned the rf frequency to many times through a region 
close to w = fi, stepping the applied magnetic field Bz 
in order to vary the Zeeman splitting 2Az. In this ex- 
periment, the molecules perform exactly as predicted by 
Eq. pSl). 



In Fig. |10| we scan a much wider frequency range, cov- 
ering the central 40 Ramsey fringes with fixed close 
to zero. Here we begin to see a departure from Eq. (23 1 



in the amplitude of the fringe pattern, which shows a 
clear decrease on either side of the centre. This hap- 
pens because the individual pulses are starting to become 
appreciably detuned, and therefore have reduced ampli- 
tude. Indeed, when the detuning reaches ±\/l5/(4r) = 
±54 kHz, the 18 /xs single pulse transition probabilities go 
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FIG. 10: The central 40 fringes of a Ramsey pattern obtained using two 7r/2-pulses of 18/iS duration, separated by T = 800 fis. 
The magnetic field is set very close to zero. Dots: Experimental points. Line: Theoretical fringes with amplitude adjusted to 
fit the data. 



to zero and the Ramsey interference has no amplitude at 
all. For a 2-level system, the envelope of the interference 
pattern is just the single-pulse lineshape, as discussed by 
Ramsey [SC. However, the spectrum of Fig. 10 involves 
the third level (u) , which makes it sensitive to very small 
magnetic fields and complicates the shape of the enve- 
lope. This detuning effect is not included in Eq. (22l, 



which assumes that the detuning is small compared with 
1/r. At present we do not have an analytical formula for 
this case, so the varying amplitude of the — cos([ri — wjT) 
curve drawn through our data in Fig. [TO]is just a smooth 
fit to the measured amplitude of the oscillations. 

The curve in Fig. [iT] shows a blow-up of the central 
Ramsey fringe (solid line), together with the fringe ob- 
tained when the high voltage leads were reversed (dashed 
line). The small shift in the phase of the fringe pat- 
tern shows that the 12.5kV/cm electric field decreased 
in magnitude by 455 ±11 mV/cm when we attempted to 
reverse it. This experiment demonstrates that Ramsey 
interferometry of the hyperfine states can provide very 
sensitive monitoring of high electric fields in a molecular 



beam apparatus. 

We turn now to interferometry using pairs of 7r-pulses. 
Whereas the 7r/2-pulses of Ramsey interferometry pre- 
pare and interrogate an rf coherence between the states 
(0) and (c), the tt pulses make and probe the superposi- 
tion of the states (-1-1) and (—1), which are almost de- 
generate. As expressed in Eq. (24l, this superposition 



evolves at the (±1) splitting frequency, leading to fringes 
of the form sin^(AzT), where hA^ = grfJ-sBz- Figure 
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shows the interferometer fringes measured in an electric 
field of lOkV/cm by scanning the applied magnetic field 
over 160 nT, with a separation of T = 800 /iS between the 
two TT-pulses. The solid line is a fit over the two cental 
fringes to the lineshape ^sin^(A2T) -t- C, where C repre- 
sents the background due to unpumped F — I molecules 
and to scattered light. Within the 2% uncertainty of the 
applied magnetic field calibration, the fringe spacing is 
found to be 44 nT, corresponding to a g- factor of gp = 1, 
as one would expect for the F — 1 hyperfine level of this 
state. The total field B^, which determines A^, is 
taken in our fit to have an adjustable offset in addition 
to the applied field. This is found to be 1.7nT in Fig. 12 
which is typical of the field that leaks through the mag- 
netic shielding from outside. While this simple function 
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FIG. 11: The central fringe of FigJTO] Solid line: fit to data 
taken with electric field as in Fig.|10[ Dashed line: fit to data 
taken with high voltage leads reversed. The evident Stark 
shift between the two indicates a change in the magnitude of 
the electric field 
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FIG. 12: Interferometer fringes obtained by scanning the ap- 
plied magnetic field, using a pair of 7r-pulses, separated in 
time by 800 /is. Dots: experimental data. Line: a fit of the 
central two fringes to the form Asin^(A2r) -|- C. 
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fits the central two fringes well, the next fringe on ei- 
ther side clearly has less amplitude. This is because the 
magnetic field detunes the two 7r-pulse transitions from 
resonance. 

Adjusting the magnetic field to reach a steep part of 
the fringe pattern, the interferometer becomes sensitive 
to any small change in the splitting 2A2 between the 
levels (±1). This can be used to monitor variations of 
the magnetic field in the apparatus. In addition, if the 
electron has a dipole moment de, then the interaction 
with the applied electric field E also contributes to A^, 
the additional amount being ±der]E/h. Here 77 is the 
enhancement factor discussed in Sec lIIBI If the electric 
field is reversed, this shift changes sign, causing the fringe 
pattern to be modulated from side to side as the electric 
field is flipped back and forth. In this way, a measure- 
ment of the signal on the side of a fringe can become very 
sensitive to the presence of a small electron EDM. 



V. ALTERNATING GRADIENT 
DECELERATION OF POLAR MOLECULES 

A. Introduction 

We have seen in section [Tv| how spectral resolution can 
be improved by increasing the time T available for coher- 
ent evolution. In a molecular beam with a velocity of a 
few hundred m/s and an interaction region 1 m long, this 
time is a few ms, corresponding to minimum lincwidths 
of several hundred Hz. For this reason, the prospect of 
decelerating and trapping molecules offers significant im- 
provements for some precision measurements, provided 
the beam flux can remain high and the inhomogeneous 
broadening due to trapping fields can be kept under con- 
trol. The basic idea of deceleration is to manipulate polar 
molecules in an electric field gradient using the force due 
to the Stark effect. After Stark deceleration was first 
demonstrated [S2], the new technology was rapidly ap- 
plied to make new measurements. Using a beam of ^'^NDa 
molecules decelerated to 52 m/s, the energies of all 22 
hyperfine levels of the (J,K)=(1,1) state were measured 
with accuracies better than lOOHz ^T. Using Stark- 
decelerated OH radicals, greatly improved measurements 
of the ground-state A-doublet microwave transitions were 
made, thus contributing to the constraint on the evolu- 
tion of the fine-structure constant over cosmological time 
[8]. Once molecules could be trapped, it became possible 
to measure directly the lifetimes of long-lived molecular 
states |53j . and even to measure the optical pumping of 
molecules by room temperature blackbody radiation |54j , 
which typically occurs on a timescale of many seconds. 

Many of the precision measurements discussed in this 
chapter make use of heavy polar molecules. Stark decel- 
eration of these heavy species is considerably more chal- 
lenging because (i) the kinetic energy to be removed is 
proportional to the molecular mass and (ii) the low-lying 
energy states are all high-field-seeking whereas the Stark 



deceleration method works best for low-field seekers. The 
first difficulty arises because the molecules formed in a 
supersonic expansion acquire the speed of the carrier gas 
into which they are seeded; thus their speed is indepen- 
dent of their mass. This difficulty could be mitigated 
by using a low-temperature effusive source such as the 
buffer-gas sources recently demonstrated |55l 156] . 

The second difficulty results from the closely packed 
rotational energy levels of a heavy molecule which causes 
all the low-lying states to be high-field seeking when 
the electric field is strong. This problem is best illus- 
trated by considering the Stark shift of a rigid rotating 
molecule of reduced mass m' , bond-length R and dipole- 
moment ^. The Hamiltonian is H = BJ^ — fl ■ E, where 
B = /{2m! E?) is the rotational constant, J is the an- 
gular momentum vector and E is the applied electric 
field. Figure [13] shows the first 16 energy eigenvalues, 
in units of i3, as a function of applied electric field, in 
units of B/fi. The field mixes states having different 
values of J but the same value of M, the projection of 
the angular momentum onto the electric field axis. Each 
energy level is labelled according to the quantum num- 
bers (J, M) that the state evolves into when the elec- 
tric field is adiabatically reduced to zero. Note that the 
states (J, M) and (J, —Ad) are degenerate for all electric 
fields, a consequence of time-reversal-symmetry. The im- 



portant point to note from Fig. 13 is that all the weak- 
field seeking states have turning points, becoming strong- 
field-seekers at high field. For example, the lowest-lying 
weak-field seeking state (1,0) has its turning point at 
an electric field of 4.9B/fi, at which point the Stark 
shift is 0.64i3. Taking YbF as an example of a molecule 
with a small rotational constant, the electric field at the 
turning point is only 18kV/cm and the Stark shift only 
0.15cm~^. This amount of energy, which can be removed 
from the molecule in a single stage of deceleration, is to 
be compared with the 682 cm^^ of kinetic energy pos- 
sessed by a YbF molecule formed in a supersonic expan- 
sion at 290m/s [461 . Clearly a very large number of stages 
would be needed to decelerate in this way. By contrast, 
the strong-field-seeking ground-state of YbF has a Stark 
shift of 10.7 cm~^ at a field of 200kV/cm, and so Stark 
deceleration in this state seems feasible. For more com- 
plex molecules, such as those of biological interest, the 
situation is even more extremely weighted in favour of 
the high-field-seekers, as discussed in more detail in |57j . 

Unlike weak-field-seekers, which are naturally focussed 
onto the axis of the Stark decelerator, strong-field-seeking 
molecules cannot be focussed using static fields. A dy- 
namic focussing scheme needs to be employed to prevent 
the molecules being pulled towards the surfaces of the 
electrodes, where the field is strongest. The alternat- 
ing gradient focussing technique solves this problem of 
transverse confinement. The molecules travel through a 
sequence of electrostatic lenses, each of which focusses 
the molecules in one of the two transverse directions and 
defocusses them in the other. The focussing and defo- 
cussing planes alternate between one lens and the next. 
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FIG. 13: The Stark shift of the low-lying energy levels of a 
rigid rotating molecule. The electric field is expressed in units 
of -B/p and the energy in units of B. States are labelled by 
the quantum numbers (J, M). 



For a subset of the molecules that enter the decelera- 
tor, namely those that lie within the transverse phase- 
space acceptance of the lens array, the net effect is to 
focus in both transverse planes. Ideally, the focussing 
and defocussing forces are linear in the off-axis displace- 
ments, and the trajectories of the accepted molecules take 
them far from the axis inside the focussing lenses, but 
close to the axis inside the defocussing lenses. Thus, the 
overall focussing is a direct result of the motion of the 
molecules, and can operate even when the defocussing 
power is stronger than the focussing. 

The first experiment to demonstrate alternating gra- 
dient deceleration of polar molecules used an array of 
12 lenses to decelerate high-field-seeking metastable CO 
molecules from 275 to 260 m/s [58]. Ground-state YbF 
molecules were decelerated from 287 to 277 m/s using a 
similar machine |59j . The transverse focussing proper- 
ties of the alternating gradient have been demonstrated 
[57j by imaging metastable CO molecules exiting from 
the decelerator. Longer machines, with more sophisti- 
cated electrode designs, should be capable of decelerating 
heavy polar molecules to rest. 



B. A model alternating gradient decelerator 

We consider a decelerator consisting of a series of elec- 
trostatic lenses whose focussing and defocussing planes 
alternate. A typical electrode structure is shown in 
a). Here, each lens is formed by applying a large 
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Fig._ 

potential difference between a pair of rods whose axes 
lie parallel to the beamline. The forward velocity of the 
molecules only changes when they pass through the fringe 
fields formed in the gap that separates one lens from the 
next. In operation, the decelerator is switched between 
three states: (i) odd lenses at high voltage, even lenses 
grounded, (ii) even lenses at high voltage, odd lenses 
grounded and (iii) all lenses grounded. We would like to 
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FIG. 14: (a) Schematic of a typical electrode structure in an 
alternating gradient decelerator. (b) Potential in which the 
molecules move, approximated using an analytical form (see 
text). Solid and dashed lines represent respectively the switch 
states (i) and (ii) of the decelerator. Positions A, B and C are 
referred to in the text, (c) At the centre of the first lens, the 
potential along x is shown by the solid line and that along y 
by the dashed line. 



work out the dynamics of molecules in such a decelera- 
tor. We therefore first generate a map of the electric field 
in the decelerator using an electrostatics solver, and then 
construct the interaction potential, using the electric field 
dependence of the Stark shift such as is plotted in Fig.[T3| 
We emphasise that this is not the electrostatic potential, 
but rather the Stark potential in which the molecules 
move. The electrode geometry shown in Fig. 14 a) is one 
of many possible geometries, some of which have been 
discussed in detail in reference j57J . The same reference 
discusses the Stark shift and shows that it varies lin- 
early with the electric field magnitude for most heavy 
molecules in the electric fields of a typical decelerator. 
The potential in which the molecules move is similar in 
form for all the electrode geometries and molecular states 
considered, though they vary in their details. Here, we 
do not specialize to a particular geometry or molecular 
state, but instead elucidate the dynamics using an ana- 
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lytical form for the potential that approximates the true 
potential produced in most cases. 

We take the potential in switch state (i) to be 



W{x, y, z) = Wo{l + b{{x/rof - {y / vof)) I {z) , (25) 
where 



ia.-ah{z' /d + L/2d) - tanh(z7d - L/2d) 
2tanh(L/2d) 



(26) 



provides a good phenomenological approximation to the 
actual z-dependence. In Eq. (25), Wq is the electric field 



at the origin, tq is a measure of the transverse aperture of 
the decelerator, and 6 measures the transverse curvature 
of the potential. In Eq. (|26l, z' = mod(z - D, 2D) - D, 



where mod(TO, n) is the remainder on division of to by n, 
D is the lens-to-lens spacing, L measures the length of a 
lens, and d measures how rapidly the potential changes 
at the exit of the lens. The potential in switch state (ii) 
is simply W{y, x,z ~ D). 

Figure 14^b,c) shows the potential for the case where 
L = 2/31?, d = l/lbD and b = 0.15. We shall use these 
parameters throughout. Part (b) is a plot of the potential 
along the beamline, W^(0,0,z), for both switch state (i) 
(solid line) and switch state (ii) (dashed line). Part (c) of 
the figure shows the transverse potentials in the centre 
of a lens, l^(a;,0,0) (solid line) and 1^(0, ?/,0) (dotted 
line). In our ideal lens, these are everywhere equal and 
opposite. 



C. Axial motion 

Consider molecules travelling down the axis of the de- 
celerator. To calculate when we should switch the poten- 
tials we introduce the concept of a synchronous molecule, 
which enters the decelerator with speed uq- We design 
the switching sequence so that this molecule is always 
at the same relative position in the periodic array every 
time the field is turned on (e.g. position A in Fig. 14(b)) 
and every time the field is turned off (e.g. position B 
in Fig. 14(b)). We refer to these fixed positions as Zon 
and Zoff- Since this molecule is always climbing potential 
hills, it decelerates as it moves along the beamline, and 
the time intervals between successive switches must be 
chosen to increase in correspondence. The required time 
sequence can be constructed using a simple algorithm. 
Consider the energy conservation equation that applies 
between the turn-on and turn-off points, Zon and Zoff , 



\Mul_^ + T^(2on) = \M{dzldtf + W{z). (27) 

Here, z is the position of the molecule at time t. Re- 
arranging and then integrating gives us the relationship 




Time (D/Uq) 



FIG. 15: The speed of a synchronous (thick red hne) and 
a non-synchronous molecule (thin blue line) as a function of 
time, using the model potential and the turn-on and turn-off 
points A and B indicated in Fig.|14[b). The inset reveals the 
step-like structure of the deceleration. 



between the n turn-off time, ioff.m Etnd the n turn-on 
time, ion,n 



n — ^^011,71 



dz 



+ 2(M/(Zon) - W{z))/M 



(28a) 

where m„ is the speed of the synchronous molecule im- 
mediately after the n*^ turn-off. The {n + 1)*'' turn-on 
is now found using 



- ^Ul_^ + 2{W{Zon) - VF(Zoff))/M, 



on,n+l 



D - (Zoff - Zon) 



(28b) 



(28c) 



Since wc now have the potential seen by the molecules 
as a function of both position and time, we can solve 
the equation of motion numerically for any molecule, 
synchronous or not. Let us define N to be the mini- 
mum number of stages required to stop the synchronous 
molecule of mass M and initial spee d ttp , through the 
relation ll^ol = Mul/{2N). 



15 



shows the re- 



Figure 

suit of such a calculation for the case where L = 2/3Z?, 
d — 1/15D and N — 80, and the turn-on and turn- 
off points are 'A' and 'B' in Fig.[l4];b). The thick red 
line shows the speed of the synchronous molecule ver- 
sus time, while the thin blue line shows how the speed 
changes with time for a molecule that has the same initial 
speed as the synchronous molecule but starts out ahead 
by 15/15. The deceleration of the synchronous molecule 
appears to be uniform, though when magnified, as in the 
inset, is actually seen to be a series of small steps refiect- 
ing the shape of the potential. These steps can also be 
seen in the main figure when the speed is low. The speed 
of the non-synchronous molecule oscillates around that 
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of the synchronous one and also consists of many small 
steps when examined in detail. If we take snapshots of 
the position and velocity of the molecules at each of the 
many turn-off times, we would no longer be able to see 
the fine structure of the motion. Between one snapshot 
and the next, the change in kinetic energy of the syn- 
chronous molecule is /S.K = MmAu = MDu/S.u/ /S.Zs = 
W{zon) — W{zos), where we have introduced the quan- 
tity to represent the change in position of the syn- 
chronous molecule between one snapshot and the next, 
and have used the fact that Az^ = D always. Since 
we have discarded the information about the fine struc- 
ture of the motion, the deceleration of the synchronous 
molecule appears to be constant, and we can convert to 
continuous variables. The above equation then becomes 
MDudu/dzs = MDd'^zJde = W{z^^ - VF(zoff)- Ap- 
plying a similar reasoning to any other molecule, de- 
scribed by the position and velocity coordinates z and 
-y, we obtain MDd^zjdt^ = W^C^on + -5) - W^z^s + z), 
where z — z — z^, and we have made the additional ap- 
proximation that w — M ^ u so that the distance moved 
by the general molecule between switching times is also 
very close to D. Subtracting the equation for the gen- 
eral molecule from that for the synchronous molecule we 
obtain an equation of motion for the relative coordinate, 
allowing us to define an effective force, i^off^ 



df^ 



M^(Zon + 5) - VK(Zon) - M^(Zoff + 2) + VK(Zoff ) 



MD 



(29) 



Introducing the relative velocity v — v — u — dz/dt, 
the left hand side of the above equation can be written 
as vdv/dz. Integrating, we then obtain 
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l/2Mi^ + (t/(z) - y(0)) = Eq, 
where Eq is a constant and 



ViS) = -JfcS d~z 



(30) 



(31) 



is an effective potential for the relative motion between 
non-synchronous and synchronous molecules. 

Figure 16 'a) shows the effective potential obtained 
from Eqs. (|25[), (26 1 and (311. The sohd line corresponds 
to Zon = — f/lOjZoff — D/3, while the dashed line has 
Zoff = ll-D/30 corresponding to greater deceleration, but 
a shallower effective potential. Since the effective po- 
tential is confining, non-synchronous molecules with in- 
sufficient energy to reach the top of the potential must 
oscillate about the synchronous molecule. We can use 
the effective force to solve for the relative motion of 
non-synchronous molecules. Parts (b) and (c) of Fig. 16 
show trajectories in phase-space obtained in this way, for 
the solid and dashed effective potentials shown in part 
(a). In each case, the thicker line separates bounded 



FIG. 16: (a) The effective potential, Eq. (31|, as a function 
of relative position. The solid line corresponds to moderate 
deceleration (turn off position 'B'), while the shallower poten- 
tial shown by the dashed line is for strong deceleration (turn 
off position 'C'). The trajectories in phase-space calculated 
using these effective potentials are shown in (b) and (c) for 
the deep and shallow potentials respectively. The thick outer 
boundaries in these plots are the separatrices. 



and unbounded motion and is called the separatrix. All 
molecules inside the separatrix will remain close to the 
synchronous molecule throughout the deceleration pro- 
cess, and the area bounded by the separatrix is the axial 
phase-space acceptance. Comparing parts (b) and (c) of 
the figure, we see that, as expected from the shallower 
potential, the acceptance is smaller when the decelera- 
tion is greater. We note that the same phase-space plots 
can be generated without making use of the effective po- 
tential, by numerically integrating the complete equation 
of motion as was done in generating Fig. |15| The trajec- 
tories in phase space then acquire the detailed structure 
shown in the inset of that figure, but are otherwise found 
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to be identical to those obtained (much more rapidly) 
from the effective potential. 

For small-amplitude oscillations about the syn- 
chronous molecule, the motion is harmonic. Expanding 



the right hand side of Eq. ( 29 ) in a Taylor series about 
z = 0, gives (Pi/dt^ - [W'izZ) - W'{zofi))z/{MD) = 
where W'{d) = dW/dz evaluated at z = a. The angu- 
lar oscillation frequency for small-amplitude axial oscil- 
lations is therefore 



W'izos) - VF'(zo„) 



MD 



(32) 



For our model potential, the solid line in Fig. 16 'a), 
the frequency can be conveniently expressed as uJz/2tt = 
0.309 uo/{y/ND). For example, if D = 30mm,and N = 
80 when uq — 300 m/s, we find (jJz/2tt = 345 Hz for turn- 
off at position 'B'. 



D. Transverse motion 

Our model potential is harmonic in the two transverse 
directions, with a curvature that varies with z. As a 
molecule travels through a lens, the curvature is very 
nearly constant until it reaches the fringe-field of the lens 
where the curvature drops rapidly to zero. To simplify 
the analysis, we make the approximation that the curva- 
ture has the constant value, W^h/r^ over the lens length 
L, is zero in the drift space of length S = D — L, and 
changes abruptly between these values. As the deceler- 
ator structure is periodic in z, it is natural to write the 
equation of motion with independent variable z rather 
than time. For a molecule with forward speed u, the 
equation of motion is 



d^x/dz"^ + K^Q{z)x = 0, 



(33) 



where Q{z) — 1 inside a focussing lens, -1 in a defocussing 
lens and in a drift region, and the spatial frequency is 



D{1) inside a defocussing lens and 0{l) in a drift region, 
with 



F{1)^ 
D{1) = 
0{l)^ 



cos(kZ) sin(K^) 
— sin(KZ) cos(kZ) 

cosh(K^) sinh(K/) 
sinh(K/) cosh(K^) 

1 kI 
1 



(36a) 
(36b) 
(36c) 



If, in moving from left to right, the molecule first trav- 
els through a region described by the matrix Mi , followed 



by a region described by M2, the matrix in Eq. (35 1 is 
simply the product , M = M2 ■ Mi . In this way, one com- 
plete unit of our alternating gradient array is described 
by M = F{L) ■ 0{S) ■ D{L) ■ 0{S), and a sequence of 
N such units is described (in a more compact but obvi- 
ous notation) by (FODO)^ . It can be shown that the 
molecular trajectories are stable if the well known condi- 
tion -2 < Tr{FODO) < 2 is satisfied (e.g. jSP]). 

We would like to know whether a molecule that enters 
the array reaches the exit. For a long decelerator, the 
above stability condition is a necessary but not sufficient 
one, since a molecule can be on a stable trajectory that 
takes it so far from the axis that it crashes into one of the 
electrodes. Rather than constructing trajectories piece- 
wise using Eqs. (36), we follow the approach first used 



by Courant and Snyder in the context of the alternat- 
ing gradient synchrotron [61' . We look for a solution to 
Eq. ( 33 ) of the general form 



x{z) = y/ei(3{z) cos{ip{z) + 5^) 

= Ai y/piz) cos ipiz) + A2 \//3(z) sin V(z) (37) 

where /3 is a z-dependent amplitude function that has 
the same periodicity as the AG array, -0 is a z-dependent 
phase, and Ci, Si, Ai and A2 are defined by the initial 
conditions. By substitution into Eq. (33l we find that 
Eq. ( 37 ) is a valid solution provided 



/ 2Wob 
Mu^rl 



(34) 



The angular frequency of the transverse oscillation inside 
a focussing lens is independent of the beam velocity u and 
is related to nhy U, — ku. 

In moving through a region of length I, from an initial 
axial position zg, the transverse position and velocity co- 
ordinates of a molecule change according to 



x/ro 



M{zo + l\za) 



za+l 



x/ro 



(35) 



^0 



Here, the dimensionless transfer matrix denoted by 
M{zo + 1\zq) takes the values F{1) inside a focussing lens. 



/3(^') 



dz' 



and 



\p" + lf3P" + ^'Q{z)f3^ = n^ 



(38) 



(39) 



To find /?, we need to make the connection to the piece- 
wise solution already given in Eqs. (36 1. From Eq. (37 1 
we have 



AlK A2K , 

x (z) — — ^(— acos Tp — sm %p) H ^(— asm tp -\- cos ip) 

VP VP 

(40) 
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FIG. 17: (a) Dashed lines are typical trajectories through a 
section of the alternating gradient array. Solid lines denote 
the envelope that bounds the trajectories of all transmitted 
molecules. The positions in the array of the focussing (F), 
defocussing (D), and drift regions (O) are indicated by the 
dotted lines. The figure shows that the beam envelope is 
largest at the centre of the focussing lens, and smallest at the 
centre of the defocussing lens, (b) Evolution of the phase- 
space ellipse through a unit of the array. 



where 



2k 



(41) 



Using Eqs. (37 1 and (40) we find the relationship be- 



tween the coordinates x,Vx/^ at position z and those at 
position z + /cell, ^coii = 2_D being the periodicity of the 
array. We make use of the periodicity constraint on (3, 
f3{z + ^coii) = and thus obtain: 



M(z + Zcoii|2) = 



cos <I> -I- a sin $ /3 sin $ 
— 7sin<i> cos $ — a sin (f> 

(42) 

where 7 = (1 -I- o?')/ (3 and $ = V'(^cgIi) is the phase ad- 
vance per unit cell. Since the integral in Eq. (38 1 is taken 
over a full period, $ is independent of z. The matrix in 
Eq. (42 1 is known as the Courant-Snyder matrix. We can 



now equate this matrix to the explicit form for the trans- 
fer matrix for one lattice unit, and so obtain (3{z). For 
example, at a distance z beyond the start of a focussing 



lens, M{z + Uik) = F{z).0{S).D{L).0{S).F{L - z). 
We then obtain <I> using the relation cos<i> = Tr{M)/2, 
and then find /3 by equating the upper right hand element 
of M to that of the Courant-Snyder matrix. 

Figure [l7|a) shows a few trajectories (dashed lines) 
calculated using Eq. (37l with kL = 1, kS = 0.5, and 
arbitrarily chosen values for and Si. The motion is a 
product of two periodic functions, one of wavelength /cell 
and the other of longer wavelength 27r/ceii/*I'- It is often 
the case that $ <C 27r, in which case the modulation with 
wavelength l^eii has a small amplitude and is called the 
micromotion, while the much longer wavelength motion 
is called the macromotion. For the case shown in the 
figure, $ = 0.387r and the separation into a micromotion 
and a macromotion is evident. If we consider a large col- 
lection of molecules, all having different values of Si and 
ei, the only constraint being that every \ei\ < e, then all 
the trajectories will be bounded by the envelope ±\//3i. 
The bold lines in the figure show this envelope. The fig- 
ure shows that the beam size modulates with the period 
of the array, reaching its maximum size in the centre of 
every focussing lens, and its minimum size in the centre 
of every defocussing lens. Since the confining and decon- 
fining forces are linear in the off-axis displacements, the 
defocussing lenses have less effect on the beam than the 
focussing lenses, this being the key to the stability of the 
alternating gradient array. As the power of the lenses 
increases, so too does the depth of modulation of the en- 
velope until, at the stability boundary where $ = tt, the 
beam size becomes zero at the centre of the defocussing 
lens. 

Using the first line of Eq. ( [37^ and its derivative to form 
the quantity + {ax + Pv^JTl)^, we find the invariant 



7^^ -I- 2axvx/n + I3{vx/nf 



(43) 



This equation defines an ellipse in the phase-space whose 
coordinates are x and Vx/^- The coordinates of all 
molecules with the same value of but different values 
of 5i lie on this ellipse. Replacing by e gives us the 
ellipse that bounds the entire set of molecules in the col- 
lection discussed above. The shape of this ellipse changes 
periodically with z, but its area is a constant, 7re. Figure 
[l7|b) shows the phase-space ellipse at various positions 
in the array. The beam is diverging as it enters the fo- 
cussing lens. Inside this lens the ellipse rotates, reaching 
its maximal spatial extent at the lens centre where its 
principal axes are parallel to the coordinate axes. The 
ellipse continues to rotate so that it is converging at the 
exit of the focussing lens, and still converging as it enters 
the defocussing lens. At the centre of the defocussing 
lens, the spatial extent is minimized and the principal 
axes of the ellipse are again along the coordinate axes. 
The beam then starts to diverge again. 

Molecules will only be transmitted if their trajecto- 
ries do not take them outside the natural boundaries 
formed by the electrodes. The characteristic size-scale 
in the transverse direction is tq and so to calculate the 
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FIG. 18: Phase-space acceptance in one transverse direction, 
as a function of kL and kS. The acceptance is shown in units 
of Qr'i. 



transverse phase space acceptance, we assume that the 
electrodes impose an aperture of size 2ro in each of the 
transverse directions. The beam characterized by e will 
be transmitted if \/]3e < tq everywhere. In particular, 
the condition must be satisfied at the point where P has 
its maximum value, /3maxi which we already know is at 
the centre of every focussing lens. The phase space ac- 
ceptance in {x,Vx) space is thus found to be Trrgfl / /3ma,x- 



Figure 18 is a density plot of the phase-space accep- 
tance in cither transverse direction as a function of the 
two dimensionless parameters that define the array, kL 
and kS. The region of highest acceptance is found near 
kL ~ 1,5 <C L, and the maximum value is 0.7447-0 
obtained at kL = 1.254, kS = 0. 

The requirement that kL ~ 1 for high acceptance 
constrains the aspect ratio of the lenses, i.e. the ratio 
L/vq. Using Eq. (34) and setting kL = 1, we obtain 
L/tq — \J~Njh. Typical values for a decelerator with a 
maximum field of 200kV/cm arc iV = 80 and h = 0.15, 
giving an aspect ratio of £/ro = 23.1. Since k is in- 
versely proportional to u, it will increase as the molecules 
slow down. To preserve the transverse acceptance, the 
lenses can be made progressively shorter so as to main- 
tain kL ~ 1. Alternatively, the alternating gradient array 
can have the structure (FO)"(£'0)", with the value of n 
decreasing down the beamline. 



E. Beyond the ideal model 

In the previous section, wc have used an idealized po- 
tential because it allows us to understand the main as- 
pects of the dynamics in a straightforward way. This 
theoretical model misses some largely undesirable effects 
that are present in reality. In the transverse directions. 



nonlinear forces are necessarily present and these re- 
duce the transverse acceptance [57l [62l ES]. The linear 
part of the force changes sign between one lens and the 
next, leading to the dynamical stability discussed above, 
but the leading order non-linear terms in the transverse 
force do not change sign between the focussing and defo- 
cussing lenses and so tend to upset the dynamical stabil- 
ity. Even when small compared to the linear terms, the 
non-linear terms can significantly reduce the transverse 
acceptance. Calculations for some typical electrode ge- 
ometries are presented in |57] . In our idealized model, 
the axial and transverse potentials are completely decou- 
pled. This cannot be achieved in any real decelerator 
because axial gradients of the electric field, needed for 
deceleration, change the dependence of the electric field 
on the transverse coordinates. In particular, the fringe 
fields at the ends of the lenses tend to increase the force 
constant in the defocussing direction relative to the one in 
the focussing direction, leading to further beam loss |57j . 
A realistic simulation of the transmission of molecules 
through an alternating gradient decelerator should use 
a three-dimensional map of the electric field magnitude 
produced by the electrode geometry of the machine, and 
the full electric-field dependence of the Stark shift. Sim- 
ulations of this kind show that the true molecular trajec- 
tories in phase space are similar to those calculated with 
our simplified model, though the acceptance volume in 
phase space may be considerably smaller. 



VI. CONCLUDING REMARKS 

In this chapter we have summarised how precise mea- 
surements on molecules are able to address important 
questions about the constancy of physical laws and the 
structure of fundamental interactions. This is a relatively 
new direction, which has emerged from a growing ability 
to prepare and manipulate molecules in pure quantum 
states. We have discussed the preparation of pulsed su- 
personic molecular beams and have explored how coher- 
ent control of the hyperfine levels can provide exquisite 
sensitivity to electric and magnetic fields. These same 
methods also provide an opportunity to search within 
molecules for more interesting effects, such as symmetry 
violations or variation of fundamental constants, which 
can be related to new physics on high energy scales. 
Deceleration and trapping are important for improving 
these experiments because they can increase the time 
available for coherent interaction with molecules from 
milli-seconds to seconds. Although we have focussed on 
decelerating heavy polar molecules using strong electric 
field seekers, the Stark deceleration method can also be 
applied to weak-field-seeking polar molecules and to Ry- 
dberg states of atoms and molecules. Alternatively, it is 
possible to decelerate and trap molecules using optical 
dipole forces. As methods for cooling, deceleration and 
trapping molecules reach ever lower temperatures and 
higher densities, and as they are extended to heavier and 
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to more diverse species, molecules will surely play an in- fundamental physics, 
creasingly important role in testing and understanding 
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